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ABSTRACT 


The objectives of this thesis are to 1) determine the 
deflections of the horizontal and vertical tank supports 
when the ship is subjected to various loading conditions, 
and, 2) determine the stresses in the hull and support 
platforms for these loading conditions. This information, 
especially the deflections of the Support platform, is 
needed to determine the forces that will be transmitted from 
the ship through the articulated vertical support system to 
the large spherical pressure vessels that carry the cargo. 
This information will permit a more complete analysis of the 
spherical tank. 


This paper details how the beam element capability of 
ICES STRUDL was used to determine the stresses and deflec- 
tions. A quarter tank section of the ship is modeled as a 
three dimensional arrangement of beams. This model is then 
subjected to loads that simulate the loads the actual quar- 
ter tank section would experience, 


The results indicate that the maximum vertical deflec- 
tions of the support platform out of a horizontal plane are 
less than a centimeter for the loadings investigated. The 
circular hold deforms out of its circular shape in all the 
loadings as one would expect. The maximum horizontal move- 
ment of any support point on the circular support platform 
is less than 3 centimeters, Relatively hign longitudinal 
stresses are found in the main deck at the transverse 
centerline of the tank, Areas of high bending stresses are 
pointed out and should be further investigated using means 
other than a beam element model due to modeling difficulties. 
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INTRODUCTION 


The first carriage of liquified natural gas (LNG) oc- 
curred in 1958 on an experimental basis with the conversion 
of an old Ci dry cargo ship into the Methane Pioneer. Suc- 
cessful commercial ventures took place in 1964 when the 
first new ships specifically designed for the transporta- 
tion of LNG were completed and put into service. At the 
present time there are approximately 20 LNG tankers in use. 

The reason for this large growth of interest in the 
marine transportation of liquified natural gas has been the 
world wide shortages that have been caused not only by in- 
creasing general demand, but also because of the demand for 
"clean" fuels which is inspired by the ecological concern 
being expressed at the present time. Natural gas is the 
cleanest burning fuel that is available at the present time 
with the exception of hydrogen. 

Natural gas is often a by-product of oil drilling 
operations and has often been considered a waste product in 
the past, especially in areas remote from consumers. Trans- 
portation of this fuel in the gaseous state would be eco- 
nomically infeasible due to the large volume required.  How- 
ever, it was determined that if the gas could be liquified, 
the required volume per unit weight could be greatly re- 
duced. 


The principal constituent of natural gas is methane 





which has a boiling point of -258.6 degrees Fahrenheit and 

a gas to liquid volume ratio of 630:1 (Ref. 1). Thus, if 
the temperature of the natural gas is reduced below -260 
degrees Fahrenheit, the same amount of fuel could be carried 
in 1/630 of the volume of the gas at room temperature. 

These low temperatures caused several technical prob- 
lems of which the principal one was probably brittle frac- 
ture in steel. Thus, the tank material must have a ductile 
to brittle transition temperature well below -260 degrees 
Fahrenheit. The materials used for the construction of the 
LNG tanks to date have been aluminum, stainless steel, Invar 
(36 per cent nickel-iron alloy) and 9 per cent nickel steel 
(Ref. 12). The problem of brittle fracture required these 
tanks to be thermally insulated from the hull of the ship 
which is constructed of mild steel. 

There have evolved two basic tank types by which LNG 
is carried—the self-supporting tank and the membrane tank. 
There are several differing designs of each basic tank with 
varying support systems, insulations, tank shapes, materi- 
als, and type of secondary barrier, if any. 

Membrane tanks are designed such that they are liquid 
tight only, with the intention that the normal loading 
forces would be transmitted through the thermal insulation 
to the hull structure on the other side. Membrane tanks 
require a secondary barrier. The secondary barrier must be 


designed to contain a catastrophic failure of the membrane 
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for a reasonable amount of time (Ref. 9). 

There are two basic types of self-supporting tanks— 
pressure vessels and free standing tanks. Pressure vessels, 
due to their shapes, are amenable to analytical stress a- 
nalysis and are constructed in accordance with pressure 
vessel requirements. Free standing tanks are basically 
tanks, usually of rectangular or trapezoidal shape, that are 
set inside the ship's hull on supports. The advantage of 
the free standing tank over the pressure vessel is that it 
can conform more nearly to the ship's hull configuration. 
The disadvantage is that the free standing tanks require a 
secondary barrier and the pressure vessels do not. As a re- 
sult of economic considerations, some recent designs have 
been of very large pressure vessels having low vapor pres- 
sures with the differences between the two types of free 
standing tanks becoming indistinguishable. The distinction 
between free standing tanks and pressure vessels varies 
among the regulatory agencies, however, the use of a vapor 
pressure of 10 psig as a demarcation limit is common. 

The United States Coast Guard, which is charged with 
the regulaticn of the marine industry from the safety as- 
pect, has recently defined four types of self-supporting 
tanks (Ref. 16). 

1. independent Pressure Vessel Tanks — These are 

pressure vessels with the vapor pressure P great- 


er than Dt20 where D is the diameter in feet and 
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where P is at least 40 psig. 

2, Semi-Independent Pressure Vessels — This category 
includes pressure vessel tanks where the vapor 
pressure is less than for the independent pressure 
vessels. 

3. Independent Structurally Indeterminate Self-Sup- 
porting Cravity Tanks — These are free standing 
tanks in which the maximum vapor pressure is less 
than 10 psig. 

4, Independent Structurally Semi-Determinate Self- 
Supporting Gravity Tanks — These are large tanks 
with vapor pressures under 10 psig that are usually 
designed using finite element and/or fracture 
mechanic techniques to pressure vessel standards. 
This class of tank generally will be required to 
have only a partial liquid tight, splash tight 
secondary barrier. 

The proposed tanks for the Technigaz design that is 


being analyzed fall into category 4. 





DESCRIPTION OF THE SHIP AND CARGO TANKS 
The Ship 
This analysis was undertaken as part of the design 
process for the Technigaz 125,000 cubic meter methane car- 
rier design. The principle characteristics of the design 


are as follows: 


Length overall 290.6 meters 
Length between perpendiculars 275.0 meters 
Breadth, molded 44.0 meters 
Depth(at main deck) 26.0 meters 
Draft 10.97 meters 
Block coefficient 0.775 


The cargo tanks consist of five spherical tanks ver- 
tically supported at their equator (Fig. 1 and 2). The 
spaces for these spherical tanks occupy a large portion 
of the deck area necessitating the requirment of very heavy. 
plating (45 millimeters thick) at the deck and sheer strake 
junction to provide longitudinal material in tension, stiff- 
ness against lateral loads from the sea, and resistance to 
buckling in the sagging condition. 

Longitudinal bulkheads, located 8.9 meters inboard 
from the outer hull except in way of tank holds, provide the 
lateral subdivision for wing tanks which are used for bal- 


last (Fig. 3). The double bottom is composed of longitu- 
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MIDSHIP SECTION - FRAME 228 
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dinal plates spaced 3.28 meters apart. In addition the 
inner and outer bottoms are strengthened by longitudinal 
stiffeners spaced .82 meters apart. 

Transversely the hull is stiffened by deep web frames 
in the double bottom and sides. The frame spacing is .85 
meters with a web frame every 2.55 meters. There are six 


decks in addition to the double bottom and the main deck. 


The Cargo Tanks 

The cargo tanks consist of five spherical tanks of 
semi-independent pressure vessel type. Tank 1 has an in- 
Side diameter of 31.6 meters and a capacity of 16,520 cubic 
meters. Tanks 2, 3, 4, and 5 have an inside diameter of 
37.5 meters and a volume of 27,610 cubic meters. This re- 
sults in a maximum capacity of 126,965 cubic meters. The 
useful capacity, using a maximum filling ratio of 98%, is 
124,450 cubic meters. 

The spherical tanks are supported, in the vertical 
direction, at the equator by a patented system of articu- 
lated parallelograms of rods and arms (Fig. 4). This sys- 
tem is attached to brackets at the tank equator and then 
extend down to the support platform. Horizontal restraint 
to pitch and roll is provided by keys and keyways. These 
keys are installed on two circles parallel to the equator 
at the approximate positions of the Tropic of Cancer and 


Tropic of Capricorn. 





E 


PIN JOINTS 


gpherica 
Tank 
ps ሩ 
; . this bar 
this bar 14 is normally 
Normally IN IN Com vession 
TENSION : 





PIN JOINTS 


Support p lat form 


FIGURE 4 
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The purpose of the vertical support system of articu- 
lated parallelograms is to allow the cargo and tank weight 
to be supported without a bending moment being introduced 
in the shell, it also allows for the thermal contraction 
of the sphere while carrying the cold cargo. The keys 
change the horizontal forces due to pitch and roll to tan- 
gential forces on the tank shell. The purpose of these e- 
laborate support systems is to allow only forces that can 
be structurally analyzed to be imparted to the spherical 
tanks. This allows the tank plating thickness to be re- 
duced to the minimum possible thickness and also allows a 
possible reduction in the secondary barrier. 

The present proposed design scheme is to fit insulation 
around the tank inside the inner hull with the inner face of 
the insulation acting as a spray shield in case of tank 
leakage, A drip pan is fitted in the bottom of the hold to 
collect any drippings and to act as a local secondary bar- 
rier. 

One of the purposes of this thesis is to determine the 
deflections of the vertical support platform and the key- 
ways in order that the forces imparted to the sphere may be 
determined prior to a structural analysis of the shpere, 
which is required before a reduction in the secondary bar- 


rier will be permitted by the regulatory agencies. 
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MATHEMATICAL MODEL 


General Discussion 

The solution method that appeared particularly suited 
to the task of solving for the deflections of the horizon- 
tal and vertical tank supports as well as solving for the 
stress levels in the various structural members was some 
type of beam or plate element solution. Recent advances 
and knowledge in this field have been extensive. Predicted 
stress levels obtained by this type of analysis have shown 
a good correlation with actual results obtained from strain 
gages on ships at sea. 

It is fortunate that these elemental solutions have 
been developed since the traditional naval architectural 
methods of stress calculations would probably have been 
inadequate for the unconventional geometry of this vessel. 
The fact that the section modulus varies with length (large 
portions of the decks are cut out to make way for the spher- 
ical tanks) makes standard calculations difficult if not 
impossible. 

There are several element type computer programs avail- 
able that have either been developed specifically for ships 
or that are general purpose. Some of these, such as ICES 
STRUDL-11, developed by Massachusetts Institute of Technol- 
ogy for use on IBM 360/370 (Ref. 7 and 8) and DAISY (Ref. 


10), developed by the University of Arizona in conjunction 
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with American Bureau of Shipping for use on both CDC 6000 
series and for the UNIVAC 1108 computers, have user oriented 
languages, Other programs exist, especially in the aero- 
space field, where the stiffness matrix is assembled by 
hand and a computer is used for the solution. 

The decision was made to use the M.I.T. developed ICES 
STRUDL-II (Integrated Civil Engineering System-Structural 
Design Language) program which is a series of computer pro- 
grams that can be used in conjunction with each other to 
solve problems in structural engineering. The STRUDL pro- 
gram was selected primarily because of time consideration. 
The time required to develop a finite element program par- 
ticularly suited to this ship structure, would have been 
prohibitive. Another advantage of STRUDL is that the lan- 
guage is easily understandable to an engineer in that the 
words and phrases in the input are such that their meaning 
is usually self-evident. 

Several different types of analysis procedures are 
currently available in STRUDL for solving framed structures 
and continuum mechanics problems. The member stiffness 
matrix in the frame analysis is computed from beam theory, 
while continuous mechanic problems are solved with the fi- 
nite element capability and the element stiffness matrix 
is computed from energy considerations. STRUDL allows for 
the mixing of members and elements provided they have the 


same number of degrees of freedom per joint and these free- 





ee 


doms correspond in type and direction. 

Analyses of ship structures have been carried out using 
beam elements as well as finite elements. Stiansen and 
Elbaloute (Ref, 15) did a thorough study of SL-7 container 
ship design using DAISY. Fenton (Ref. 4) used the finite 
element capability of STRUDL to obtain stresses in a cata- 
maran cross structure. Zoller (Ref, 17) combined beams and 
finite elements in a tanker study with satisfactory results. 

It is very important to use the analysis procedure that 
is best suited for the problem since the number of elements 
or members required for even a coarse node network of a 
complex structure can become extremely large with the stiff- 
ness matrix computer core requirements exceeding the capac- 
ity of the computer. 

The decision was made to use beam elements in this 
study. The main reason being that the complexity of the 
structure, caused by the circular hold for the spherical 
tanks, required a fairly fine mesh with all joints having 
six degrees of freedom to get any meaningful deflections. 
This would have required an extremely large number of ele- 
ments if the finite element capability was utilized, since 
plate bending elements would have to be superimposed on 
plane stress elements to account for loads in all three co- 
ordinates. Any combination space frame and finite element 
would be completely out of the question for the same reason. 


The stiffness analysis carried out by STRUDL is a 
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linear, elastic, static, small displacement analysis. The 
procedure used requires the specification of the model ge- 
ometry, the member and joint loads, the member properties 
and various applicable constants, such as, Young's modu- 
lus and Poisson's ratio. The analysis used treats the joint 
displacements as unknowns. When a space frame is specified 
as the structural type, as it was in this problem, there 
are six unknowns at each joint—three displacements and 
three rotations. 
The stiffness analysis is comprised of the following 
steps: 
1. Consistency checks of members 
. Generation of the stiffness matrix- 
ㆍ Processing of member loads 


c Stiffness matrix assembling 


. Solution of the matrix 


2 

3 

4 

5. Processing of the joints 

6 

7. Joint displacement processing 
8 


. Processing of member stresses 
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Beam Model 

The area of the midship section around Tank 3 was the 
section selected for study. This area was chosen because 
the maximum bending moment occurred in this region for the 
loading conditions that were to be analyzed. This area 
would then give the maximum amount of deflection of the 
support platform as well as the largest stresses—the two 
parameters that were to be determined. 

In order to keep the bandwidth of the stiffness matrix 
as small as possible and still get valid results it is ad- 
vantageous to limit the size of the model. First it was 
possible to limit the model size by the use of symmetry of 
both the structure and the loading of the ship about the 
centerline of the vessel for the cases being investigated. 
This cut the size of the model in half. Furthermore, with 
approximate symmetry of the moment curve and internal struc- 
ture of the ship about frame 228, it was decided that suf- 
ficient symmetry did exist to further reduce the model. 
Thus the model could be reduced to a quarter tank section 
extending from frame 195 to frame 228 with the necessary 
boundary conditions at frame 195 to simulate the remainder 
of the vessel. 

All joints in the transverse plane of symmetry through 
the center of the tank and in the longitudinal plane of 
symmetry along the centerline of the vessel must be sup- 


ported. STRUDL assumes that all joints specified as sup- 
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port joints to be rigidly supported according to structural 
type with no displacements or rotation unless specifically 
released in a certain direction. Thus, to insure an accu- 
rate representation of the actual ship, it was necessary to 
release the constraints in the Y and Z direction as well as 
the rotation constraint about the X axis for joints in the 
transverse plane of symmetry. In the longitudinal plane of 
symmetry, it was necessary to release the force constraints 
in the X and Z direction in addition to the rotational con- 
straint about the Y axis. The origin, the center of the 
tank hold in the double bottom, remained completely fixed 
to provide a reference point. 

The structural members of the ship were idealized, as 
indicated previously, by beams. The object of this modeling 
process was to accurately represent the ship, while at the 
same time trying to keep the number of beam elements to a 
minimum because of computer space considerations. The model 
that was finally settled upon idealized two actual ship 
decks as one deck in the model. In the horizontal plane 
the flanges of the beam element were the deck plating and 
the vertical members formed the web of the beams. The 
effective breadth of plating was calculated for bending in 
both the transverse and longitudinal directions. Examples 
of beam elements formed from the double bottom and support 
platform are shown in figures 5 and 6, 


The vertical beam elements were handled in basically 
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DERIVATION OF BEAM ELEMENT CROSS-SECTIONS 
FROM CELLULAR DOUBLE BOTTOM 
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SEE FIGURE 2 
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FIGURE 6 


DERIVATION OF BEAM ELEMENT CROSS SECTIONS 
FROM SUPPORT PLATFORM DECKS 
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the same way, Generally the beam elements consisted of tee 
bars welded to plating. This was true for the outer verti- 
cal hull, the bulkhead at the inner web frames, and the cir- 
cular support platform. Figure 7 depicts the geometric 
shape of the beam model. Figures 8 through 16 show the in- 
dividual joints and member numbers for the various joints 
and bulkheads of the model. 

The major problem that developed in the modeling pro- 
cess involved the circular support platform. One charac- 
teristic of beam elements is that the ends must always be 
ata joint. Thus it was necessary that the vertical mem- 
bers at least butt on a longitudinal or transverse member 
of the double bottom. These joints have numbers 75, 115, 
156, 198, 251, and 253 in the model(Fig. 8). This restric- 
tion prohibited the vertical members from being equally 
spaced as in the actual ship. The equal spacing in the 
ship design was accomplished by the use of vertical stiff- 
ner plates in the double bottom under the vertical members. 

The model did not take into account such details as 
the rounded biige or rounded deck edge, however, the effect 
of these idealizations were judged to be negligible when 
compared to mesh used. 

It should be pointed out here that the point of con- 
nection of the beam elements is at the neutral axis of the 
element. This is the reason for the small difference in 


the actual beam and depth when compared to the model scale 
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FIGURE 7 


THREE DIMENSIONAL VIEW OF THE BEAM MODEL 
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FIGURE 8 


BOTTOM DECK OF BEAM MODEL 
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FIGURE 9 


SECOND DECK OF BEAM MODEL 
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THIRD DECK OF BEAM MODEL 
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FIGURE 11 
FOURTH DECK OF BEAM MODEL 
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FIGURE 12 
MAIN DECK OF BEAM MODEL 
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FIGURE 13 
OUTER SKIN OF BEAM MODEL 
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INNER WEB BULKHEAD OF BEAM MODEL 
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FIGURE 15 
TRANSVERSE SECTION AT FRAME 195 
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FIGURE 16 
PROJECTION OF TANK HOLD 
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beam and depth. This point is illustrated by the following 
example. The depth of the ship is 26.0 meters. The scale 
depth of the model is 23.95 meters. This is caused by the 
fact that the neutral axis of the double bottom beams are 
1.15 meters from the bottom plating and the neutral axis of 
the beams modeling the main deck are .90 meters down from 
the main deck. 

The location of the average neutral axis was used as 
the location of all the neutral axes of the beam elements 
that would be located ina straight line on the actual ship. 
This can be best illustrated using the outer web frames 
where the thickness of the shell plating changes as the 
height is increased. This change in plate thickness causes 
the neutral axis of each beam element to shift slightly. 
The member properties of the individual members were cal- 
culated using their own neutral axis. However, in assem- 
bling the model geometry, the average neutral axis was used 
allowing the members to lie ina vertical straight line. 

There are disadvantages to using the beam model that 
was outlined in the preceeding paragraphs. One disadvantage 
is that the beam elements are connected at the joints by a 
point. The actual area of connection however may extend a 
considerable distance from the point that results from the 
intersection of the neutral axis. This is especially true 
when plates are being modeled as beams which was the case 


with this model. The other principal disadvantage of this 
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model was the result of modeling two actual decks as a 
Single deck. This would cause a vertical beam that had a 
support at each end to be modeling a member that actuaily 
was supported at three places along its length. Along the 
same line a vertical member that is subjected to a uniform 
or linear force would only be supported at the ends while 
the actual member would be supported in at least three 
places. The result of these two disadvantages would be that 
the local bending stresses in the individual members would 
be inaccurate. Figures 17 and 18 illustrate these two dis- 
advantages of the beam model. 

However, it should be pointed out that these disadvan- 
tages would only affect local bending stresses to any great 
extent. Overall results, such as joint displacements and 
longitudinal stresses, would be virtually unaffected. In 
addition the bending stresses that resulted from this model 
would give an indication of areas that should be looked at 
more closely—possibly with a fine mesh finite element 
model of that small area using the results of this beam 
model to determine the boundary conditions. 

The geometric model was entered in the STRUDL program 
by first numbering each joint and giving its coordinates in 
the global X Y Z space. The orientation of the members was 
determined by the coordinates of the joints that were its 
end points. Each member has a local coordinate system 


associated with it for the easy identification of input and 
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BEAM MODEL 
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The light lines in the top figure indicate material from 
which the member properties of the beams were calculated. 
The bending of the vertical plate about the Y axis is not 
accurately modeled by the two vertical beams since the 
junction is actually a straight line and in the model it 
consists of two points. 


FIGURE 17 
MODELING ILLUSTRATION 
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Figure (A) shows a beam supported at four points with a 
uniform load applied. Figure (B) illustrates how the same 
structure would be modeled. As can be seen by the dashed 
line showing the relative magnitude and sign of the bending 
stresses, the model gives a local bending stress that does 
not represent the actual case. However when the deflections 
of the joints are considered the modeling is accurate since 
a single beam of the model would have a strength equivalent 
to two of the actual beams. 


FIGURE 18 


MODELING ILLUSTRATION 
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output. The local positive x axis goes from the starting 
joint to the end joint. The one remaining degree of free- 
dom for the member was rotation about the principal axes. 
This last degree of freedom was specified by the angle beta. 
The reference position for the beta angle (beta = 0) is 
with the local axis z in the same direction as the global Z 
axis and with the local axis x parallel to the global Y 
axis. The method of determining the beta angle for other 
cases is outlined in Reference 7. 

There were many locations on the actual ship where 
additional stiffners were added to the shell plating between 
structural members. To account for their effect the pris- 
matic cross-sectional area of these stiffners were "blended" 
into the plate to produce an effective plate thickness 
(Sample Calculation 1). 

It was not necessary to determine the effective width 
of plating for any of the members since no buckling of any 
member was anticipated. However, it was necessary to de- 
termine the effective breadth of the various beams, both in 
the y and z directions in order that a realistic moment of 
inertia and section modulus could be calculated for each of 
the local member coordinates. 

Originally the plan was to determine the effective 
breadths using Schade's curves (Ref. 13 and 14). However, 
the use of these curves required knowledge of the loading 


conditions prior to the calculation of member properties. 
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The intention of this work was to do several loading con- 
ditions, (loaded, ballasted, hogging, sagging, etc.) and 

it would have been completely infeasible, time wise, to go 
back and recaiculate the member properties of all the beams 
for each loading condition. Instead, the standard design 
criteria for mild steel was used where the effective breadth 
of plating was equal to sixty times the plate thickness. 
This criterion results from the fact that the effective 
breadth of plating acting in association with a structural 
member, is influenced by the yield strength and the modulus 
of elasticity of the material (Ref. 11). This influence can 


be calculated by the formula 


effective breadth = 2VE/F ቲ 
where E = modulus of elasticity (psi) 


e tensile yield strength 


The above design criteria results from the fact that 2VE/F, 
equals 60 for mild steel. 

STRUDL requires that the cross-sectional axis, the 
torsional rigidity coefficient, and the moments of inertia 
about the local y and z axes be specified for the stiffness 
analysis of a space frame. Since the sectional stresses in 
the members were desired it was also necessary to specify 
the section moduli about the local y and z axes, The shear 
areas in the Y and Z directions were also supplied in order 


that the analysis would also include shear deformations. 
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A computer program was written to calculate the above 
member properties to avoid tedious and repetitive hand cal- 
culations. The program was capable of calculating the mem- 
ber section properties of prismatic beams of five or fewer 
elements. The procedures used by the computer program are 
illustrated by an example in Sample Calculation 2. 

The beam element model was constructed in such a way 
that it could be expanded to handle non-symmetrical loadings 
about the center, or in other words, the expanded model 
would be capable of handling the case when the ship is ina 
heeled position. The final quarter tank section model was a 
space frame consisting of 414 beam members and 208 joints. 

Information that was requested as output from the 
STRUDL program was: (1) the forces acting on both ends of 
the member, (2) the reactions at support joints, and (3) the 
displacements of each joint. In addition section stresses 
of the members were requested for each end and the center of 
the member. The section stresses were broken down into 
axial and bending stresses from which the maximum and mimi- 
mum normal stresses were calculated. The shear stresses at 
the three sections Were calculated by dividing the shear 


forces by the shear area, 
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THE ANALYSIS—HOGGING CONDITION 


General Discussion 

The first case to be investigated was with the ship in 
an upright position, fully loaded with the wave crest amid- 
ships, the trochoidal wave used by Technigaz to compute the 
shear and bending moment curves had a length equal to the 
ship's length-between-perpendiculars which was 275 meters. 
The wave height was equal to 0.03 of the wave length or 8.4 
meters. The depth of water at frame 228 was then calculated 
by adding the wave height to the stillwater draft. The 
Smith effect was also included in these calculations. 

An interesting load curve results when a ship carries 
cargo in large spherical tanks as with this ship. The 
spherical tank is supported at various points around its 
equator, thus the loads are transmitted to the hull of the 
ship in a circle, The load curve would then have extremely 
hign loads per unit length in the areas where the tanks are 
adjoining and, fairly low loads in the area of the trans- 
verse centerlines of the tanks. Thus, for this particular 
ship the load curve would have high spikes at frames 61, 
129, 195, 261, 327, and 389. This loading curve cannot be 
accurately approximated by a series of straight lines indi- 
cating a uniform load over a certain segment of the vessel's 
length as is the normal design practice. The load curve for 


this particular type of vessel must be either accurately 
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plotted and a graphical solution obtained or some type of 
curve-fitting computer solution performed. 

The shear and moment curves provided by Technigaz were 
calculated using the dynamic weight of the cargo. The 
dynamic loads in this region were calculated by determining 
the vertical accelerations at tank 3. The vertical acceler- 
ation included the acceleration due to gravity, heave accel- 
erations and pitch accelerations, the total of which was 
13.21 meters/second. These accelerations were determined 
from Series 60 data that gave results indicating a heave of 
+7 meters in 10.0 seconds and a pitch of 15.7 degrees in 
10.0 seconds. The dynamic force was then calculated by 
multiplying the mass of the tank and LNG times the ratio of 
the total acceleration to the acceleration of gravity. This 
ratio, which was the number of "g" forces felt at tank 3,was 
1.4. The resulting dynamic load due to the tank and cargo 
was 20,200 metric tons. 

Information supplied by Technigaz indicated that the 
steel weight, excluding the spherical tank, was 6,472 metric 
tons between frames 195 and 261, the area of interest. A 
quarter tank section then had a steel weight of 1,618 metric 
tons. 

The only additional forces that had to be applied to 
the model were the boundary forces which will simulate the 
remainder of the ship. As is indicated by Figure 19, the 


after portion of the ship may be simulated by the applica- 
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tion of a bending moment of 484,940 meter-tons and a shear 
force of 1,221 tons at frame 195, the free end. In addition 
it was necessary to apply a shear force at frame 228, the 
transverse plane of symmetry, of such sign and magnitude as 
to maintain the model in rigid body equilibrium. Or, in 
other words, it is necessary that the sum of the forces in 
the Z direction be equal to zero. An alternative method 
used by some investigators in this field to insure rigid 
body equilibrium is to fix the support joints in the trans- 
verse plane of symmetry in the Z direction. This method is 
not satisfactory in this particular case since knowledge of 
deflections of the members is the primary reason for the 
analysis and for this particular ship there are large Z 


direction deflections in the plane of symmetry. 


First Run—loads 

The beam model was to be subjected to basically four 
different types of forces. These were boundary conditions, 
loads due to water acting on the hull, loads due to the 
cargo, and loads due to the steel hull weight, The question 
was, how should the various forces be distributed to accu- 
rately represent the actual case. 

The actual steel weight of the quarter section of the 
model was distributed equally among the joints with a joint 
load of 7.78 metric tons in the negative Z direction. This 


was considered a valid distribution since there were a 
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greater number of joints in locations where there was a 
large amount of structural steel. The concept of a uniform 
force in the negative Z direction applied to each member was 
completely infeasible from the standpoint of additional pro- 
gramming and computer time required, although, this would 
have been the ideal solution, 

The dynamic load of the cargo was 20,200 metric tons of 
5,050 tons per quarter tank section. This load was applied 
as a joint load at each of the vertical members of the sup- 
port platform (third deck of the model). The load was 
equally distributed among the joints resulting in a load of 
505.0 tons at all the joints except joints 22 and 338 which 
were in the planes of symmetry and consequently had one half 
the applied load of the other joints or 252.5 tons. Sample 
Calculation 3 illustrates the method of calculation of the 
applied loads. 

The forces due to the water acting on the bottom and 
side of the hull were calculated as if the depth of the 
water remained constant over the entire quarter tank sec- 
tion. This assumption was slightly in error in that the 
depth was slightly less at the ends of the model due to the 
fact that the wave crest was near the center of the model. 
However, it was felt that this effect was negligible due to 
the trocoidal nature of the wave crest. The forces on the 
Side of the ship were idealized as linearly varying loads on 


the vertical members. The bottom forces were applied to the 
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model as uniform forces—half of the total bottom forces 
were applied to the longitudinal members and half to the 
transverse members. A ratio of ship dimension to model 
dimension was used in the calculation of these forces so 
that the force actually experienced by the ship is applied 
to the model. 

The remaining loads to be applied were the boundary 
conditions. The bending moment of 484,940 meter-mtons at 
frame 195 was simulated by coupled forces. The neutral axis 
of the model at this point was calculated to be about 10.5 
meters from the double bottom. Above the neutral axis a 
uniform load was applied in the positive X direction to the 
main deck between the inner and outer web frames. In addi- 
tion a linear load was applied to the inner and outer webd 
bulkheads from the main deck to the neutral axis, with the 
third and fourth decks carrying a uniform load equal to the 
value of linear distribution at the intersection of the webs 
and deck members. The other half of the couple was applied 
below the neutral axis in the negative X direction. This 
was accomplished by imposing a uniform load across the 
entire half width of the double bottom, with a linearly de- 
creasing load on the webs, starting at the double bottom up 
to the neutral axis. The second deck carried a uniform load 
equal to the linear load applied to the web at its joints. 
See Figure 20 for a graphical representation of the distri- 


bution. 
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FIGURE 20 
FRAME 195 BOUNDARY CONDITIONS 
HOGGING CONDITION—-RUN I 
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No forces were applied to members nearer the center- 
line than the inner web with the exception of the double 
bottom. The reason for this was that structural members 
do not contribute to the longitudinal strength of the vessel 
due to their discontinuous nature in way of the holds for 
the spherical tanks. The inner web was the structural sec- 
tion that was closest to the centerline of the vessel and 
maintained its longitudinal continuity although it did form 
part of the tank hold. 

The shear load at frame 195 was 1,221 mtons or 610.5 
mtons for the quarter tank section. This was equally dis- 
tributed as a joint load of 61.05 mtons to each of the 10 
joints that lie in the intersections of the plane of frame 
195 and the inner and outer bulkheads of the wing tanks. No 
shear force was applied to the vertical bulkheads near the 
centerline because of their large lightening holes and poor 
vertical strength when compared to the inner and outer web 
bulkheads. 

The remaining, boundary condition is the shear force at 
frame 228. This force was determined by summing the verti- 
cal forces applied to the model. The negative of this 
summation then gave the force necessary to keep the model 
in equilibrium. For this particular loading condition this 
force was a negative 1,389 mtons which was distributed as 
joint loads at nodes in the inner and outer web bulkheads 


at frame 288. See Figure 21. 
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FRAME 228 BOUNDARY CONDITIONS 


HOGGING CONDITION—RUN I 
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First Run—Results 

The first attempts to get results from this model were 
unsuccessful due to the large amount of core and facilities 
required by the model. Some of the subprograms that make up 
STRUDL have job control statements inherent in them that 
allow for a certain number of disks to be in operation when 
that subprogram is being used. This particular model re- 
quired more discs than were provided by the subprogram. 

Thus it was necessary to enter input statements which over- 
rode some of the job control statements of some of the sub- 
programs and allowed the computer to utilize sufficient 
facilities. 

Results for the first run were finally obtained and the 
deflections appeared satisfactory when looked at on a macro- 
scale. However, a closer examination of the section 
stresses indicated several problem areas. A major problem 
occurred at frame 195. The deflections of the fourth deck 
at the inner web in the positive X direction were large over 
a beam length as indicated by Figure 22. This resulted in 
all members eminating from joint 408 to have extremely high 
bending stresses. Furthermore, this large deflection in the 
X direction was transmitted to the circular hold bulkhead at 
joint 283 causing the circular support platform to elongate 
with resulting high stresses in all members in the vicinity. 

The elongation of the circular hold acting in conjunc- 


tion with the water forces on the side of the ship, caused a 
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FIGURE 22 
DEFLECTIONS OF FRAME 195 
HOGGING CONDITION—RUN I 





deflection in the negative Y direction of the circular 
support platform bulkhead at the narrow area of the ship's 
hull near the transverse centerline of the tank (frames 216 
to 228), This effect was most pronounced at the second deck 
level and caused high bending stresses in the vertical mem- 
bers in this area. 

À second major problem area was the main deck, again 
at the narrower sections between frames 219 and 228. High 
stresses occurred in the horizontal members due to the high 
axial stresses combined with the high bending stresses. 

The third major problem area occurred in the double 
bottom at the junction of the double bottom and the vertical 
strength members of the support platform again in the vicin- 
ity of frame 228, The joints at this junction had large 
negative Z defiections, which caused large bending stresses. 

A fourth problem area was in the vicinity of the origin 
in the double bottom. Members eminating from the origin had 
a positive Z deflection for a radius of about 6 meters in- 
stead of the expected negative deflection. 

Figures 23 and 24 show the longitudinal stress distri- 
bution at frames 195 and 228 respectively. Approximately 
25 per cent of the members of the model had total normal 
stresses that exceeded the yield stress of mild steel. At: 
this point a complete review of the model and all calcula- 
tions was undertaken. 


The first discovery was that an error of about 5 per 
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FIGURE 23 
LONGITUDINAL STRESSES—FRAME 195 
HOGGING CONDITION—RUN I 
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FIGURE 24 
LONGITUDINAL STRESSES—FRAME 228 
HOGGING CONDITION—-RUN I 
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cent had been made in the calculation of the bottom water 
pressure forces acting on the longitudinal members. This 
resulted in the fact that the sum of the forces in the 2 
direction was not zero. This imbalance in forces caused a 
small amount of rigid body motion which would then have 
caused the unexpected deflections in the vicinity of the 
origin outlined above as the fourth problem area, 

Conversation with representatives of American Technigaz 
brought to light the fact that the dynamic weight of steel 
had been used in the calculation of the shear and bending 
moment curves instead of the static weight. This factor had 
not been taken into account for the first run. This cor- 
rection would increase the negative Z direction joint loads 
that simulated the steel weight. A summation of the forces 
in the Z direction then indicated that the shear force re- 
quired at frame 228 to prevent rigid body motion could then 
be reduced from 1,389 mtons to 742 mtons. Consequently, 
correction of this error would help solve the problem of 
high stresses in the double bottom at the intersection of 
the double bottom and the vertical members of the circular 
support platform which was described earlier as problem 
three. 

It was obvious that the extreme elongation of the cir- 
cular hold and the inconsistent X direction deflections in 
frame 195 were caused by the boundary forces that were used 


to simulate the bending moment. Thus it was necessary to 
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determine whether the distribution of forces that was used 
was correct and if not, what was the correct distribution. 
Obvious methods, such as building a plastic scale model of 
two or three sections of the ship, loading the model and 
getting the longitudinal stress distribution at frame 195, 
were considered. Another method that could have been used 
to obtain the desired information was a large macro-mesh 
model consisting of plate elements which could be developed 
using the finite element capability of STRUDL. These two 
methods of solution had to be discarded because of time 
considerations. A literature search was then conducted to 
try to find longitudinal stress data from experiments using 
plates with circular cutouts. 

Papers were found in which experiments had been per- 
formed on plates subjected to a uniform tension force at its 
ends with a series of circular cutouts in which the cutouts 
had been reinforced with a combing. This would have approx- 
imated very closely the main deck of this ship. However, 
the papers that contained this work only presented stress 
results for a transverse section through the center of the 
circle. This is understandable since the highest stresses 
would occur at this point. No information regarding the 
stress distribution in other parts of the plate was pro- 
vided. 

The only information that was readily available was 


Kirch's solution for a central circular hole in an infi- 
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nitely wide plate subjected to a pure tension stress (Ref. 


5). The results of his solution were: 


r < 2 : 2 " 
-0 ےہ‎ © G 5 5 EO e 
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Է E r + 
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(1 + 2 > (1 + 3 1 ) cos 20 
Ի È 
where g = tension stresses to which plate is subjected 
ره‎ > radial stresses 
da > tangential stresses 
Bern radius of hole 
r = radius at which calculation is being made 
6 = angle r makes with longitudinal line through 


center of the circle 


Subsequent to the calculation of the radial and tangential 
stresses through a transverse section 20.34 meters ahead of 
the center of the circle it was possible to calculate the 


distribution of the longitudinal stresses by the equation: 


o. 一 G, ԸՕՏՑ + og Sin 


where ايه‎ longitudinal stress 


The result of the calculations utilizing the above 
equations is shown in Figure 25 with the stresses indicated 


as a fraction of the stress at what would be the edge of the 
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main deck. 

There are several important differences between the 
actual main deck plating and the idealized plate with the 
circular cutout. First, in the actual situation the main 
deck is not infinite in the horizontal plane as was the 
above solution. Secondly, Kirsch's solution did not take 
into account the fact that the actual case had an inner web 
bulkhead running longitudinally or that it had a vertical 
wall around the circumference of the hole. In addition the 
solution for the infinite plate with the single circular 
cutout did not take account of any change in the stress 
distribution that a series of holes in the longitudinal 
direction would cause. 

In spite of the many differences between the actual 
main deck and the plate with the circular cutout and with 
no way of quantitatively determining the effects of these 
changes, it was felt that the calculated distribution would 
be helpful in determining the correct boundary conditions 
that should be imposed on the beam model. 

It was readily apparent from a comparison of Figures 
23 and 25 that, in no way, did the distribution of longitu- 
dinal stresses that was the result from the first run com- 
pare with the distribution obtained from Kirsch's solution. 
For example, the longitudinal stress in the plate calcula- 
tion at the inner web bulkhead was 19 per cent of the load 


at the point corresponding to the side of the ship while the 
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results from the first run indicated that the stress at the 
inner web was 88 per cent of that at the main deck edge. 
Thus it was obvious that a more realistic application 
of the bending moment boundary conditions was necessary to 
obtain any meaningful results. This new distribution of 
forces would have to put a greater load on the part of the 
main deck that is longitudinally continuous as well as a 
greater share on the outer hull. Such a refined approach 
would tend to reduce the problems encountered on the fourth 
deck at the intersection of the circular hold and the inner 
web which resulted in an extreme elongation of the circular 
hold. A more realistic distribution would also tend to re- 
duce the high stresses that occurred at the narrow portion 


of the main deck at frame 228. 


Second Run—Loads 

The decision was made to run the STRUDL program a sec- 
ond time for the hogging condition making the changes to the 
loads that were necessary. The first correction made was to 
change the joint loads due to the ships steel weight from 
7.78 metric tons/joint to 10.89 metric tons/joint. This new 
weight was the result of using a "g" force of 1.4 as was 
previously discussed. This change also made necessary a 
change in the shear force at frame 228. A summation of 
forces in the Z direction indicated that the shear force per 


joint should be reduced to 74.2 mtons/joint as indicated in 
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Figure 26. The force due to the water acting on the long- 
itudinal members of the double bottom was also corrected for 
this run. 

The largest problem for the second run was to determine 
the distribution of longitudinal forces that would accurate- 
ly simulate the bending moment and give a reasonable distri- 
bution of longitudinal stresses in the main deck. After a 
careful reanalysis it was decided that no forces should be 
applied to the horizontal members making up the second, 
third, or fourth deck. The reason for this was that these 
decks had almost no longitudinal strength due to the large 
lightening holes in them. The elimination of forces on 
these members would be a conservative estimate when the ship 
is looked at as a whole, since the decks would assume a 
small portion of the longitudinal load. In order to try to 
obtain a longitudinal stress distribution in the main deck 
a linearly varying load was applied to the transverse mem- 
bers of the main deck of frame 195 between the inner and 
outer web frames. The value of the force at the inner web 
was half of that at the edge of the ship as indicated in 
Figure 27. In addition, a linear load was applied to the 
vertical member of the inner and outer web bulkheads, with 
the maximum value being the force on the main deck or the 


double bottom and a zero force at the neutral axis. 


Second Run—Results 


The results from the second run seemed much more rea- 
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sonable than had those of the first run. The positive Z 
direction deflection of the double bottom near the origin 
had completely disappeared with the elimination of the rigid 
body motion that had been caused by the miscalculation of 
the longitudinal forces. However, as can be seen from 
Figure 28 the longitudinal stress distribution on the main 
deck at frame 195 still did not come close to the distribu- 
tion obtained for the exact solution of the infinite plate 
with the circular hole. The results from Run II indicate 
that the stress at the inner web is approximately 54 per 
cent of the longitudinal stress at the edge of the main 
deck. This 54 per cent represented an improvement over the 
88 per cent obtained in the first run but still falls short 
of the desired result of 19 per cent obtained from the in- 
finite plate solution. It was obvious that the forces 
applied to the main deck that are part of the couple used to 
simulate the bending had to be redistributed again, 

The problem on the fourth deck where the inner web 
joined the circular hold wall was greatly improved, however, 
the stress in some members was still excessive. The bending 
stresses remained the reason for this problem with members 
427, 428, and 429 being pulled out of the circular shape by 
member 409. The poor distribution of the bending moment 
couple forces that was discussed in the previous paragraph 
would be one of the primary reasons. 


As can be seen from Figure 29 the longitudinal stresses 
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in this run still exceeded the yield stress in the main deck 
at frame 228, This problem was discussed at a conference 
with Technisaz representatives. At this conference it was 
discovered that there had been a communication problem con- 
cerning a revision to the blueprints from which this work 
was done. The revision concerned an additional longitudinal 
bulkhead under the main deck that had been added for struc- 
tural purposes., Runs I and II were made using 20 milli- 
meters as the thickness of the bulkhead. The correct thick- 
ness of this bulkhead was 45 millimeters. Thus it was nec- 


essary to make a third run to obtain valid results. 


Third Run—Loads 

Most of the loads for this run remained the same as for 
the second run, the only exception being the forces making 
up the upper half of the couple that was simulating the 
bending moment at frame 195. The distribution of these 
forces was changed such that the linearly varying force be- 
tween the inner and outer web had a value at the inner web 
that was one third of the value at the deck edge. The ver- 
tical members in the inner and outer webs were also subjec- 
ted to an X direction linearly varying load that had its 
maximum value at the main deck and had a value of zero at 
the neutral axis. The boundary conditions at frame 195 are 
shown graphically in Figure 30. The boundary conditions at 


frame 228 remain the same as shown in Figure 26, 
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FIGURE 30 
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In addition, the properties of members 501 through 518 
were recalculated to account for the error in the bulkhead 


thickness. 


Third Run—Results 

The results from the third run were better than the 
first or second runs. The shape of the longitudinal stress 
distribution in the main deck at frame 195 (Fig. 31) con- 
formed more closely to that of the infinite plate solution. 
The value of the longitudinal stress at the inner web was 
45 per cent of that at the deck edge. Although this is not 
the value obtained from the infinite plate model, time con- 
siderations as well as the high cost of each run indicated 
that further experimentation with the boundary conditions 
simulating the bending moment on the upper half of the ship 
was unwarranted. This decision was based on two considera- 
tions. The first was the uncertainty of the actual longi- 
tudinal stress distribution - the infinite plate had many 
drawbacks as was previously cited. The second consideration 
was that the distribution of the forces on the main deck 
would have very little effect on the vertical deflections of 
the support platform - the primary desired result. Figures 
33 through 37 show the deflections obtained for the third 
run. 

The maximum vertical deflection of any point in re- 


lation to any other point on the circular platform for the 
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first, second, and third runs were 1.0 cm, 0.9 cm, and 0.6 
cm, respectively. Even with major changes in boundary condi- 
tions (between run I and run II) and changes in the longitu- 
dinal strength of the main deck (between run II and run III) 
this deflection difference remained fairly constant. 

Further support for the argument that the distribution 
of forces would have little effect on the vertical deflec- 
tions can be seen by looking at the magnitude of the verti- 
cal deflections between runs II and III in which the verti- 
cal loads were constant. The smallest vertical deflection 
on the support platform occurred at joint 22 and was -6.5 cm 
and -6.6 for runs I and II respectively, while the largest 
occurred at joint 338 and was -7.4 cm and -7.2 cm respec- 
tively. The conclusions being that changes in the distri- 
bution of forces simulating the bending moment and even 
changes in the longitudinal strength of the main deck have 
little effect on the vertical deflections of the support 
platform. 

The same is, naturally, not true for the transverse or 
longitudinal deflections since the boundary conditions will 
greatly affect the longitudinal elongation and subsequent 
transverse contraction of the circular support platform when 
the ship is in a hogging situation. The maximum deflection 
in the X direction of the support joints occurred in the 
longitudinal direction and was due to the elongation of the 


circular support platform. This joint deflection was 4.4 
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cm, 3.3 cm, and 2.5 cm for the first, second, and third 
runs respectively. 

The largest Y direction deflection occurred at the 
transverse centerline of the tank. This contraction of the 
circular platform in the transverse direction was directly 
linked to the elongation in the longitudinal direction. The 
maximum magnitude of this transverse deflection was 4.1 cm, 
2.8 cm, and 2.1 cm for the three runs. 

The third run had the most realsitic longitudinal dis- 
tribution. Thus, the maximum expected deflection of any 
joint normal to the circular shape for this loading condition 
was 2,5 cm, It was felt that this number was conservative 
and could be used for further work on the preliminary de- 
sign. 

The areas of high stresses that occurred in the third 
run were in basically the same areas as in the two preceed- 
ing runs, although the magnitude of the stresses were lower. 
The maximum longitudinal stress occurred at the narrow por- 
tion of the main deck as would be expected. The longitudi- 
nal stresses when combined with the bending stresses re- 
sulted in a total member stress that was greater than yield 
for the beams in this area. 

Vertical members in this same area near the narrow por- 
tion of the circular hold between the lower two decks had 
combined stresses that exceeded the yield stress. The 


primary cause of this was the bending stresses which result- 
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ed from a combination of loads involving the transverse con- 
traction of the circular hold area and the loads due to the 
water pressure acting on the side of the hull. In addition 
vertical members of the circular support platform at the 
ship's longitudinal centerline indicated high bending 
stresses due to the longitudinal elongation of the circular 
support platform. The same elongation produces high bending 
stresses in members 427 and 428, 

It was felt that the fact that beam elements were being 
used to model plates in the outer skin of the ship was the 
cause of the high bending stresses in this region. The 
linear force distribution simulating the water pressure 
caused the members to bow in, resulting in high stresses. 
However, in the actual case there would be another point of 
support in the middle of the beam, cutting the amount of 
bowing as was previously explained in the mathematical model 
chapter. 

However, it was the opinion of the writer that steps 
should be taken to insure that the bending stresses due to 
the elongation and contraction of the circular hold are sub- 
stantially smaller than the results indicated for this 
model. The same is true for the normal stresses in the main 
deck at the narrow portion of the hull near the transverse 
centerline of the tank. 

The increase in the area of members 501 and 510 should 


have brought about a proportional decrease in the longitudi- 
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nal stress level of the main deck at frame 228, This stress 
reduction could not be checked directly for each member be- 
cause the distribution of boundary forces in the main deck 
at frame 195 had been changed between run II and run III. 
However, the axial force on each member was one of the out- 
puts of this program. These forces were used as correction 
factors to obtain the stress level of the members if the 
same distribution had been used for run II as was used for 
Bun 111. 


This was accomplished by using: 


F = oA 
— መው و‎ 
and Իշ = و0۸‎ ԻՉ = OA, 
since A, = A, 
ت‎ vi 
F,/0, - 2 
ժշ Հ 0, (F3/F2 ) 
where J, = actual longitudinal stress obtained for run Il 


kk 


Jd, = longitudinal stress that would have been ob- 
tained for run II if the boundary conditions 
of run III had been used 


Բշ = axial force in member for run III 
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F. = axial force in member for run III 


հշ = cross sectional area of member for run II 


Since the stresses and area can be considered inversely 
proportional if a constant force is applied, the graph of 
the stress versus area of a member would be approximately 
linear. Thus, the information from runs II and III can be 
used to determine the area required for any desired stress 
level when the vessel is in the hogging condition. See 


Figure 38. 
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THE ANALYSIS—SAGGING CONDITION 


General Discussion 

The last case to be investigated was with the ship in 
an upright position, fully loaded with the wave crests at 
the end of the ship and the wave hollow at amidships. Again 
the wave height was taken to be 0.03 of the wave length or 
8.4 meters. The draft at frame 228 was then calculated by 
subtracting the wave height from the stillwater draft, this 
resulted in a draft of 6,80 meters, 

The shear and moment curve for this loading condition 
were also porvided by Technigaz. The bending moment at 
frame 195 was -309,800 meter-tons while the shear load was 
296 mtons as indicated by Figure 39. Since the dynamic 
loads would be slightly less than the static loads for the 
sagging condition, the static loads were used to give con- 


servative results. 


First Run—Loads 

Again the beam model was subjected to boundary loads, 
loads due to water acting on the hull, cargo loads and loads 
due to the steel weight. The static steel weight of 7.78 
mtons/joint was applied to each joint of the model. 

The cargo load was 14,428 metric tons or 3,607 metric 
tons per quarter tank section. The load was equally distri- 
buted among the joints resulting in a load of 360.7 metric 


tons at all the joints except joints 22 and 338 which were 
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in the planes of symmetry and had one half the applied load 
of the cuter joints or 180.4 metric tons. 

The forces due to the water acting on the bottom and 
side hull were calculated as if the depth of water remained 
constant over the entire quarter tank section. The forces 
on the side of the ship were idealized as linearly varying 
loads on the vertical members. As in the hogging case the 
bottom forces were applied to the model as uniform forces— 
half of the total forces to the longitudinal members and 
half to the transverse. See Sample Calculation 3. 

The bending moment was simulated by coupled forces 
taking the neutral axis to be 10.5 meters from the double 
bottom. The bending moment of 309,800 meter-mtons was dis- 
tributed as it was for Run III of the hogging case. Above 
the neutral axis a linear load was applied in the negative 
to the transverse member of the main deck between the inner 
and outer frames. The magnitude of the force at the inner 
web was one third of the force at the deck edge. A linear 
load was also applied to the inner and outer bulkheads from 
the main deck to the neutral axis. The other half of the 
couple was applied below the neutral axis in the positive X 
direction. A uniform load was imposed across the entire 
half width of the double bottom, with a linearly decreasing 
load on the webs, starting at the double bottom up to the 
neutral axis. See Figure 40 for a graphical representation. 


The shear load of 296 mtons for the whole cross section 
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or 148 mtons for the quarter tank model was equally distri- 
buted as a joint load of 14.8 mtons acting in the positive Z 
direction on each of the ten joints that lie in the inter- 
section of frame 195 and the inner and outer web bulkheads. 
The one remaining boundary condition, the shear at frame 
228, was determined by summing the vertical forces applied 
to the model and then changing the sign of this summation 
thus causing the model to be in equilibrium in the vertical 
direction. This shear at frame 228 was determined to be 
1598 mtons. This resulted in a joint load of 159.8 mtons 
for each of the joints in the inner and outer web bulkheads 


of frame 228 (Fig. 41). 


First Run—Results 

Generally, the deflections and stresses for the sagging 
case were less than for the hogging case. As would be ex- 
pected the directions of the longitudinal and transverse 
deflections were the opposite of those in the hogging case. 
The circular tank hold contracted in the longitudinal direc- 
tion and elongated in the transverse direction. 

The shape of the longitudinal stress distribution in 
the main deck conformed almost exactly to the distribution 
obtained in the third run of the hogging case although the 
sign was different because the deck was in compression for 
this case (Fig. 42). The longitudinal stress at the inner 


web was 44 per cent of the stress at the deck edge. This 
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FIGURE 42 
LONGITUDINAL STRESSES—FRAME 195 
SAGGING CASE—RUN I 
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LONGITUDINAL STRESSES—FRAME 228 


SAGGING CASE—RUN I 
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result is very close to the 45 per cent obtained for the 
hogging case. If the main deck had been an infinite plate 
subjected to pure compression or tension the distributions 
would have been exactly the same. 

The deflections in the sagging case were of a smaller 
magnitude than the other runs as can be seen from Figures 
L4 through 48. Figure 46 indicates that the maximum deflec- 
tion of any point in relation to any other point on the cir- 
cular support platform was 0.3 cm. The maximum deflection 
of the joints in the support platform carrying the cargo 
load was 0.7 cm and occurred at the transverse tank center- 
line due to the transverse elongation. The maximum contrac- 
tion of the support platform occurred near the ship's longi- 
tudinal centerline and had a magnitude of 0.5 cm. 

There were only four members that had a total normal 
stress ir excess of the yield stress. Two of these members, 
609 and 627, are vertical members in frame 195 and the high 
local bending stresses are due to the bowing caused by the 
application of the boundary conditions. The high stresses 
are completely unrealistic and are of no concern in the 
actual ship. The other two members, 672 and 687, have high 
stresses due to bending. This bending was caused by longi- 
tudinal contraction of the circular tank hold in the case of 
member 672 and the transverse elongation in the case of mem- 
ber 687. The stresses in these two members bears further 


investigation. However, it was the opinion of the writer 
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FIGURE 45 
SECOND DECK DEFLECTIONS 
SAGGING CASE-—-RUN I 
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FOURTH DECK DEFLECTIONS 
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MAIN DECK DEFLECTIONS 
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that these high bending stresses are the result of the beam 


model and are not realistic. 
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CONCLUSIONS 


1, The ICES STRUDL program can be effectively utilized 
in the preliminary design phase to determine stresses and 
deflections of ships of unusual form for which traditional 


design criteria is unsatisfactory. 


2, The beam element model used gave good results for 
the deflections of the circular support platform in the 
vertical direction. These results were basically indepen- 
dent of end effects due to the application of loads simu- 
lating the remainder of the ship. The maximum vertical de- 
flection of any point in relation to any other point on the 
circular support platform was 0.6 centimeters for the con- 


ditions investigated. 


j, The results for the magnitude of the transverse and 
longitudinal deflections of the circular platform were in- 
fluenced by end conditions, particularly the distribution of 
forces simulating the bending moment. Upper limits were 
obtained for these magnitudes using a conservative distri- 


bution of forces. 


4, The longitudinal stresses were also influenced by 
the end condition. However, it was possible to obtain an 
upper limit on the maximum value using a conservative dis- 


tribution of boundary conditions. This maximum longitudinal 
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stress occurred in the main deck at the narrow portion of 
the hull as would be expected. The stress level in this 
area was of sufficient magnitude to warrant consideration of 


the use of high strength steel in this area. 


5. The local member bending stress results obtained 
did not 2 reflect the actual case. This is due to 
the inability of beam elements to accurately model plates 
on a local scale. However, indications of areas with high 
bending stresses were obtained and these should be subjected 
to further investigation. Possible reduction of cutout sizes 
or additional flange reinforcement at the holes are possible 


cures if needed. 
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RECOMMENDATIONS 


This thesis investigated the stresses and deflections 
for the hogging and sagging condition with the ship in an 
upright mode. It is recommended that additional work be 
carried out to determine the ships reaction to rolling and 
to torsion. The torsional stresses incurred when this ship 
is heading obliquely to a wave would be appreciable due to 
the small amount of main deck area, 

A further recommendation is that the end effects be 
eliminated in any future work with this model or any other 
single tank section model. This can be accomplished by any 
of three ways: 

(1) The model could be expanded in the longitudinal 
direction, thus eliminating the end effects on the quarter 
tank section of interest. 

(2) A scale, strain guage equipped model could be con- 
structed so that the correct distribution of longitudinal 
stresses in the main deck and web bulkheads at frame 195 can 
be determined. 

(3) A macro-scale finite element program could be 
developed that would give as output the longitudinal stress 


distribution at the main deck. 
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Sample Calculations 
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SAMPLE CALCULATION #1 
EFFECTIVE THICKNESS OF PLATE 


Double Bottom - Transverse Cross Section 


T,= /6 mm 


v 
+ Ty? 18,5 mw 





ا و د పలాల‏ 


stiffner spacing = 0.82 meters 
frame spacing = 3.28 meters 
cross sectional area of 400 X 14 HP stiffners = 6125 mm^ 


cross sectional area of 400 X 15 HP stiffners = 8450 mm^ 


At = area of stiffner y (no, stiffner between frames) 
frame spacing 


t = 16 + 9122 X 3 = 21.6 mm 
l effective 3280 
0 
Ե = 18,5 + X 3 > 26.2 mm 
effective 3200 
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SAMPLE CALCULATION #2 


MEMBER PROPERTIES 


Double Bottom - Transverse Cross Section 





ا س 
m T, 2 26,2 mw‏ 328 1,3 
:4« 2 


local coordinates 


The double bottom longitudinals may be idealized as 


beams as indicated by the solid lines above. 


AX = cross-sectional area 


= (3.28) (.0216) + (2.50) (.0185) + (3.28) (.0262) 
- 1 


AY = shear area in Y direction = area of flanges of beam 
= (3,28) (.0216) + (3.28) (.0262) 
- 7ه‎ m< 


AZ = shear area in Z direction = area of web of beam 
一 (2, 50) (40185) 
= 6 2ا0‎ m< 
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IX = torsional rigidity 
3 


(1/3) (3.28)(.0216)3 (3.28)(.0262)% (2.50)(.0185)? 


000036 m 


The effective breadth of the plating is equal to 60 times 
the thickness of the plate 


in the y direction 


Ly = 60(.0216) = 1.296 m 
eff 
Lo = 60(.0262) = 1.572 m 
eff 
in the z direction 
L4 = 60(.0185) = 1.11 m 
eff 


It is necessary to determine the y and z neutral axes in 
order to calculate the moments of inertia and section mod- 
uli for the y and z axes. 


y neutral axis - it is necessary to use I, and 
eff 


Lo to calculate the y neutral axis 
eff 





YNA 


YNA 


TY 


o 
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21.6 mm 


YN A 





*—— 1. 2 4599. —— À 
es 


= distance from bottom of member to neutral axis 


distance from the y neutral axis to the extreme section 


= J (areas) (moment arms 
> areas 


— (1.296)(,0216)(2,50) 4 (2,50)(.018 
— (1.296)(.0216) 4 (2,50)(.0185) 4 (1. 52271. 6265) 


1,148 m 


1 


moment of inertia about y axis (using L, and Lj,  ) 
eti ㆍ eff 


(1/12) (.0185) (2.50)?  (.0216)(2.50 - 1.107)? 
+ (,0262) (1.572) (1.107) 


21 ኃሪሪ m 


section modulus about y axis 


۱ 
© 
NO 
~J 
00 
= 
ኒልን 
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Z = distance from the z neutral axis to the extreme section 


The z neutral axis is a vertical line through the horizontal 
centers of the three sections of the member. 


TZ = moment of inertia about Z axis (use L4 ) 
eff 


3 
= (1/12) (.0216) (3.28) + (1/12) (.0262) (3.28) 


= 6 m^ 


IZ 
7 
02 

.6 


= ,0857 ոշ 


Տե = 


H 
: 
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SAMPLE CALCULATION #3 
The following load calculations are examples using 


values taken from Run III of the hogging case. 
NOMENCLATURE AND VALUES 


I draft 
= 15.2 meters 


P = density of water 
= ء‎ 6 mtons/m' 


M = moment at frame 195 
= 484,900 mtons-meter 


Q = shear at frame 195 
= 610 mtons 


W = static steel weight of quarter tank section 
= 74,55 tons/meter 


= = acceleration of gravity 
Eo m/sec“ 
a, = acceleration due to heave 
= 2.76 n/sec* 
a, = acceleration due to pitching 
= ,04 m/sec” 
n = number of joints in beam model 
= 208 joints 
c = static load of cargo (quarter tank section) 
= 3607 mtons 
r, = ratio of actual ship's bean to model's beam 
በ 2፡00 
u et 


o 052 


Բլ 5 ratio of actual quarter tank section length to model's 
length 


0۵+ تی 





EM 


Forces Due to Water Acting on the Hull 


Ք = pressure on bottom hull = (T) )۾(‎ 
ا‎ O20) (sy 2) 
15.6 mtons/em* 


| 
5 
The water pressure loads were applied to the vertical 
members as indicated below, 


— — VATER 


15,2 ንጻ = 


(P)(member width) 
au 6) (2,55) 
39,78 mtons/cm 


= 
ዘ 


9 for the vertical members at frame 228 is half of the 


above value due to symmetry conditions. 


The loads on the bottom were divided such that half of the 
loads were applied to longitudinal members and half were 
applied to the transverse members. 


F = (P)(member width) (r,) 


w 
eo) er 55) (1.1059) 
21.06 mtons/m 


վ 


F, for transverse members 1, 10, 19, 28, 38, 49, and 60 is 


one half the above value due to symmetry. 


The forces on the longitudinal member were calculated in a 
similar fashion using 


በ (P)(member width) (r,) 
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Loads Due To Cargo 


un un ግ 5 LD nan es z ም " La = = E e y^ Ese 
ከወ 65... 1105!ሮ OF Torce felt by the quarter tank due to 


Are mom NUS ec‏ صا رج 
mur M 020,3‏ 


- 
| 


dynamic cargo load = (c)(g) 


co 


= (3607) (104) 
= 5050 mtons 


There are 11 support joints, two of which only carry half of 


౧ ౧౯ ine others, due to symmetry, Thus: 


E ,5 4 .5) 


if 
ایا‎ 
ርጋ 
LA 


mtons/support joint 
Loads Due To Steel Weight 


- (wig) 
۷ 


) 
= (24,55)(21,7)(1.4) 


= e265 mons 


cr 
RI 


w YW 
2265/ 208 
10,89 mtons/ joint 


/ 


1 
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BOUNDARY LOADS - FRAME 195 


The horizontal neutral axis was calculated to be located 
10.5 meters above the bottom deck. The loads simulating 
the bending mement were assumed to be distributed as inci- 
cated in Figure 30. The moment that had to be applied to 
the quarter tank model was one half of the actual moment of 
the whole ship at frame 195. The forces applied above the 
neutral axis would form half of the couple. Thus to deter- 


mine FL and ES 





21/7 couple 
1 F 2F L F.L 
Wi t L 1 L 2 2 
F L 
+ 2 2 
where 7 = distance between the inner and outer webs 
= 6,27 m 
L, = distance between the neutral axis and the 


main deck 
13.45 m 


1 


բ բ P. 
SEES ENDS 2(8.37) (13.4) * (8.37) (13.45) + ”ፓ ա 


F 
+ 313.45)" 


mtons/meter‏ 779.94 رہ 
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ԽՆ 
M _ b 2 
LE Fpl T3) (Ly) + (2) 5 (3 Ly) 


where եշ = width of double bottom 
జ 
Ly er distance from double bottom to neutralMaxis 


10.5 
484,940 _ 2 
5 b O 2 (lO e ስ 6270) 
FE 445.6 mtons/meter 
The shear load at frame 195 for the quarter tank section 


was 610.5 mtons. This was distributed equally between the 


ten joints shown in Figure 30. Thus 


Q/10 


= աա = 61.05 mtons 


ርጋ 
| 


Boundary Loads - Frame 228 


In order to calculate the shear forces at frame 228 to pre- 
vent rigid body motion of the model it is necessary to make 
the summation of the forces in the vertical direction equal 


to zero. 


F summation of Z direction forces 


Z 


[ 


Q -= EP 一 ት + water on bottom hull 


610.5 - 2265 - 5050 + 7437 
742 mtons 





1 
HA 

I E 
=š 


| 
~J 
= 


۰٥0677۴ 
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APPENDIX B 


Computer Input for Run III 
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